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Abstract 



The NA44 collaboration has measured charged kaon and pion distributions at 
midrapidity in sulphur and proton collisions with nuclear targets at 200 and 
450 GeV/ c per nucleon, respectively. The inverse slopes of kaons are larger 
than those of pions. The difference in the inverse slopes of pions, kaons and 
protons, all measured in our spectrometer, increases with system size and is 
consistent with the buildup of collective flow for larger systems. The target 
dependence of both the yields and inverse slopes is stronger for the sulphur 
beam suggesting the increased importance of secondary rescattering for SA 
reactions. The rapidity density, dN/dy, of both K + and K~ increases more 
rapidly with system size than for tt + in a similar rapidity region. This trend 
continues with increasing centrality, and according to RQMD, it is caused 
by secondary reactions between mesons and baryons. The K~/K + ratio 
falls with increasing system size but more slowly than the p/p ratio. The 
7r~/-7r+ ratio is close to unity for all systems. From pBe to SPb the K + /p 



data suggest that as larger nuclei collide, the resulting system has a larger 
transverse expansion, baryon density and an increasing fraction of strange 
quarks. 

PACS numbers: 25.75.-q 13.85.-t 25.40.ve 
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I. INTRODUCTION 



Nucleus-nucleus collisions at ultrarelativistic energies create hadronic matter at high 
energy density. Enhanced production of particles containing a strange quark, such as kaons, 
may indicate formation of a state of matter in which the quarks and gluons are deconfined 
]I|fJJ^,|]. This is because while the production of s-quarks is generally suppressed in hadronic 
collisions, they are light enough to be abundantly produced at temperatures above the 
deconfmement phase transition, and rapidly reach chemical equilibrium via gluon-gluon 
fusion (see |Q for a review). Alternatively there may be a duality of quark matter and 
"conventional" hadronic explanations for the strangeness enhancement observed in recent 
data || . The K~ and K + yields together provide a sensitive probe of the spacetime evolution 
of heavy-ion reactions. Since K~s have a large annihilation cross-section with neutrons, their 
yield is sensitive to the baryon density. The K~ and K + distributions may also hint at the 
degree of thermalization achieved; and their transverse mass spectra allow detailed study 
of rescattering and collective expansion effects. The pions are less sensitive to collective 
expansion because of their small mass. However because they are the most numerous of 
the produced particles they give information on the total entropy produced in the collision. 
We present charged kaon and pion measurements using the NA44 spectrometer from pBe 
collisions (to approximate pp), p-nucleus and nucleus- nucleus interactions. This allows a 
systematic study as a function of the size of the central region and different conditions in 
the surrounding hadronic matter. 

The distributions of kaons at midrapidity provide a sensitive probe of the collision dy- 
namics and constrain the assumptions of event generators. We compare our data to the 
RQMD model, version 1.08 RQMD is a microscopic phase space approach to modeling 
relativistic heavy ion collisions based on resonance and string excitation and fragmentation 
with subsequent hadronic collisions. RQMD also contains "ropes" formed by overlapping 
strings, whose sources are color octet charge states 0. These increase the number of strange 
antibaryons and antiprotons but have little effect on kaon yields or slopes. 

II. EXPERIMENT 

The NA44 experiment is shown in Fig. [l|. Three conventional dipole magnets (Dl, D2, 
and D3) and three superconducting quadrupoles (Ql, Q2, and Q3) analyze the momentum 
and create a magnified image of the target in the spectrometer. The magnets focus particles 
from the target onto the first hodoscope (HI) such that the horizontal position along the 
hodoscope gives the total momentum. Two other hodoscopes (H2, H3) measure the angle of 
the track. The hodoscopes also measure the time-of-flight with a resolution of approximately 
100 ps. Particle identification relies primarily upon the third hodoscope. 

The momentum acceptance is ± 20% of the nominal momentum setting. The angular 
coverage, with respect to the beam, is approximately ± 5 mrad vertically and -5 to +78 
mrad in the horizontal plane when the spectrometer is at 44mrad with respect to the beam. 
When the spectrometer is rotated to 131mrad, the horizontal coverage is 77 to 165 mrad. 
Only particles of a fixed charge sign are detected in a given spectrometer setting. Four 
settings are used to cover the midrapidity region in the pr range to 1.6 GeV/c. Fig. |2| 
shows the acceptance of the spectrometer for pions and kaons in the y — Pt plane for the 
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4 and 8GeV/ c momentum settings when the spectrometer axis is at 44 and 131 mrad with 
respect to the beam. 

Two Cherenkov counters differentiate kaons, protons and pions and reject electrons. An 
appropriate combination of CI and C2 was used for each spectrometer setting to trigger 
on pion or kaon/proton events. Particles are identified by their mass-squared, constructed 
from the track momentum and time-of- flight, in addition to the Cherenkov information. 
This is shown in Fig. |3] for the 4GeV/ c setting. Pions are clearly visible after kaons and 
protons have been rejected by requiring a signal in the first Cherenkov counter. Kaons and 
protons are clearly separated after pions have been rejected using the same Cherenkov. After 
mass-squared and Cherenkov cuts the residual contamination of the data is less than 3%. 

A scintillator (To) is used to trigger on central events in sulphur-nucleus collisions by 
requiring a large pulse height (high charged particle multiplicity). The pseudorapidity cov- 
erage of T is roughly 1.3 to 3.5. For proton beams, T provides the interaction trigger 
by requiring that at least one charged particle hit the scintillator. More details about the 



spectrometer are available in [10,11 



III. DATA ANALYSIS 

The data samples after particle identification and quality cuts are shown in Table |. Also 
shown for each data set is the target thickness and the centrality, expressed as a fraction 
of the total inelastic cross-section. In order to construct the invariant cross-section, the 
raw distributions are corrected using a Monte Carlo simulation of the detector response. 
Simulated tracks are passed through the full analysis software chain and used to correct 
the data for geometrical acceptance, reconstruction efficiency, particle decay in flight and 
momentum resolution. Particles are generated uniformly in rapidity and according to an 



exponential distribution in transverse mass, m T = ^Jp^ + m 2 , with the coefficient of the 
exponent determined iteratively from the data. 

The absolute normalization of each spectrum is calculated using the number of beam 
particles, the target thickness, the fraction of interactions satisfying the trigger, and the 
measured live time of the data-acquisition system. For the SA data, the centrality selection 
is determined by comparing the pulse height distribution in the T counter for central and 
unbiased beam triggers. For pA systems, the fraction of inelastic collisions producing at 
least one hit in the interaction (To) counter is modeled with the event generators RQMD 
and FRITIOF fl2]l3|| . The resulting centrality fractions are indicated in Table Q. 



The proton-nucleus data are corrected for non-target background. The largest corrections 
are 15%, 12% and 12% for K + , K~ and 7r + respectively, for pBe collisions. This correction 
does not affect the shape of the distribution, only the absolute cross-section. No correction 
is needed for the central nucleus-nucleus data. The cross-sections are also corrected for the 
particle identification and trigger inefficiencies. 

The invariant cross-sections, measured in the NA44 acceptance, are generally exponen- 
tials in transverse mass. This allows us to characterize the distributions by two numbers: 
the inverse slope, and the total particle yield. A detailed description of the estimated sys- 
tematic errors on particle yields and inverse slopes is given in [TT| . These were derived from 
checks on the momentum scale, acceptance corrections, trigger and lifetime efficiencies, see 
Table 0. 
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IV. RESULTS 



The invariant cross-sections for kaons and pions from pBe, pS, pPb, SS, and SPb colli- 
sions are shown in Figs. f| and [| as a function of rrij- — m. The transverse mass distributions 
are generally described by exponentials in the region of the NA44 acceptance, that is:- 

1 Ed3(T = Ce ~(m T -m)/T^ /jx 

a d 3 p 

where C is a constant and T is the inverse logarithmic slope. The pion spectra are not 
exponential for all rnj- however as discussed in Section [V|. 

The inverse slope parameters obtained by fitting the spectra to Eq. [I] are given in Table [D.1 
and plotted in Figs. |6| and [5]. The inverse slope parameters for both K + and K~ increase 
with system size. If the SPb sample is split into two centrality bins, 11 — 6% and 6 — 0%, the 
inverse slopes for both K + and K~ show no centrality dependence. The inverse slopes of 
kaons from RQMD are smaller than the data. The shape of the RQMD spectra for kaons is 
concave in the NA44 acceptance, resulting in slopes that are smaller than the data. Further 
study is necessary to understand why this is so. The inverse slopes of pions show a weaker 
system dependence than kaons but increase with %. 

The inverse slopes of our ir + , K + and proton |ll]] spectra are shown in Fig. |8] versus 
system. The inverse slopes of pions increase only slowly with system size, while those of 
kaons and protons increase more rapidly. The difference in slopes increases with system size. 
This is consistent with transverse flow, which should be stronger for larger systems. 

The rapidity density, dN/dy, was calculated by integration of the normalized rriT — m 
distributions, corrected as described above. The fitted slope was used to extrapolate to 
high ttiti beyond the region of measurement. The statistical error on this extrapolation was 
calculated using the full error matrix from the fit of Eq. |T] to the ttit spectrum. 

The K + and K~ dN/dy values from data and RQMD are plotted in Fig. ^| and the 
data (including the pion dN/dy values) are given in Table |V|. The uncertainties in dN/dy 
introduced by the uncertainty in the slope of the my spectrum, the Cerenkov trigger effi- 
ciency and the combining of the spectra from the two settings are included in the errors. 
The number of kaons found at midrapidity increases rapidly with the size of the colliding 
nuclei. Comparing (pPb/pS) with (SPb/SS) shows that the target dependence is stronger 
in nucleus-nucleus collisions than in p-nucleus collisions. 

Fig. |H] shows the ratio of dN/dy for K~ and K + for the various target-projectile systems 
studied. Note that the systematic errors listed in Table [TJ cancel in this ratio. The ratio for 



pBe agrees with K~ / ' K + from pp collisions measured at the ISR The ratio drops 

slowly with system size but this is not explained by the lower yfs of the SA data. Only a 3% 
change in K~ / K + was observed for pp collisions when changing y/s from 31 to 23GeV/c 



14j| . The ratio also decreases for more central collisions. Note the it /tt + ratio is 1.0 ± .15 



for all systems. RQMD reproduces the K~ / ' K + ratio. 

Fig. |ll| shows the yield of positive pions and kaons and their ratio for y m 3 vs system. 
For pions, like kaons, the target dependence is stronger in nucleus-nucleus collisions than in 
p-nucleus collisions. The K + /n + ratio is constant for pA collisions but increases by about 
70% in central SA interactions. The ratio also increases with centrality. 

A similar effect has been observed by the E802 collaboration for SiAu collisions at 
14.6GeV/c, where the K + /tt + ratio is a factor of 3 ± 1 larger than for pp collisions of a 



5 



similar energy fT6| . RQMD was able to give a reasonable explanation of this result JT? 



in 



terms of secondary production of K + by meson-baryon interactions, see also ||. A similar 
result was obtained with ART |Tj|, while in ARC [^] the K + enhancement was caused by 
A A interactions. RQMD predicts quite well the measured features of the dN/dy of pions 
and kaons, but slightly overestimates the K + /ir + ratio in pA collisions. 

For SS collisions, isospin conservation implies that the total yield of Kg over the full 
phase space should equal the average of the K + and K~ yields. (A similar argument holds 
for pions). Comparing our yields with those of NA35 p0| , we find excellent agreement for 
kaons and reasonable agreement for pions [pT{] . However the NA35 kaon inverse slopes are 
significantly larger 



22 



The NA36 and WA85 collaborations have measured the inverse 
slopes of charged kaons from SPb and SW collisions respectively and we are in good agree- 

Although WA85 does not give absolute yields, their K~ / K + 

2.3-2.9. 



23.24 



ment with their results 

ratio for y = 2.3 — 3.0 from SW |24j] is very close to our value for SPb at y 



V. DISCUSSION 



Fig. [7] shows that the tt + inverse slopes increase with increasing p?- Cronin et al [25 
parametrized their pA pion spectra in terms of pp spectra as follows 

—(pA) = A^— (pp). (2) 

This worked well for many different nuclei. It was found that a depended only slowly on 
rapidity but strongly on p?, running from 0.7 at low p? to greater than 1 for very high p^. 
The low p T value of alpha was close to the value 2/3 expected for an opaque nucleus while 
a = 1 implies full participation of the nucleus. In [^6| a > 1 was understood in terms of 



parton scattering. The HELIOS collaboration extended Eq. ^| to nucleus-nucleus collisions 
|27| by writing 



^ {AS)={ AB^^ M . (3) 

Using Eq. ^ to compare our SPb and pBe ir + data, we find that a = 0.8 at pt = and rises 
to 1.0 by about 1.2 GeV/c. This is in agreement with recent results on tt° production by 
WA80 ||28|| . HELIOS observed slightly smaller values of a when comparing negative hadron 



production from SW and pW . This may be due to their more backward rapidity range 
where the density of particles is lower. 

A prerequisite for creating a quark gluon plasma is to produce a hot dense nuclear system 
in heavy ion collisions. Such a system would be likely to undergo transverse expansion which 
would affect the inverse slopes of the particle spectra. The effect of such a velocity boost 
would be more pronounced on heavier particles. Fig. |] shows that the difference in inverse 
slopes between pions, kaons and protons increases with system size. This suggests that 
the effect of transverse flow becomes stronger for larger systems. NA44 has previously 
reported an increase with system size in the inverse slopes of kaons, protons and antiprotons 
produced at mid-rapidity in symmetric collision systems ||29[| . These data extend this trend 
to asymmetric systems. Since the effect of the boost increases with mass, the inverse slopes 
of pions do not increase much with system size. 
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The target dependence of dN/dy is stronger in nucleus- nucleus collisions than in p- 
nucleus collisions probably because additional particles are made in secondary collisions of 
the produced particles. The increase of the K + /n + ratio with system size suggests that the 
fraction of strange quarks is increasing. In RQMD, rescattering increases kaon production by 
50% for SS, mainly because of associated production of K and A in meson-baryon collisions 
@. However, the total multiplicity remains essentially the same if rescattering is turned on 
or off in the code. 

Fig. [12| shows the K~ / K + ratio versus our p/p ratio [|ll|] for various systems at y=2.3-2.9. 



If we have a system of quarks in chemical and thermal equilibrium, then ratios of particles 
and antiparticles can be described by ratios of the chemical potential of the constituent 
quarks to the temperature f3(| . Since the quark content of the proton is uud and K + = us 
(K~ = us), then K~ / K + = e 2fls ^ T (p/p) 1 ^ 3 where /i s is the chemical potential of strange 
quarks and T is the temperature. If /i s /T changes slowly from pBe to SPb then we would 
expect K~ / K + oc (p/p) 1 / 3 . The data are close to this simple form. As the colliding system 
gets larger, the baryon density increases reducing the fraction of u quarks and driving down 
the K~/K+ ratio. 

Fig. [13| shows the ratios K + jp , K~ jp and their geometric mean \J K+ p ^ versus sys- 



tem. The K + /p ratio decreases while K~ jp increases and y K+ p .p stays constant. Again 
assuming chemical and thermal equilibrium, this implies that the number of strange quarks 
is proportional to the product of the numbers of u and d quarks. Such behavior is consistent 
with the increase of the K + /tt + ratio shown in Fig. |TT[ This suggests that in SA collisions 
strange quarks are still suppressed somewhat with respect to u and d quarks. In thermal 
language, the chemical potential of the u and d quarks is not yet equal to the mass of the 
strange quark. 



VI. CONCLUSIONS 

Several features of the data suggest a strong build up of rescattering and hence energy 
density as we collide larger ions. These are: 

• The stronger target dependence of yields and slopes for SA collisions than pA collisions; 

• The increase of the pion inverse slopes with m-r; 

• The increase in the K + /tt + ratio with system size. 

If particles undergo many rescatterings, pressure builds up and this causes radial flow. The 
increase of the inverse slopes with particle mass and the increase of the difference in slopes 
for larger systems supports the existence of radial flow. However the extraction of a flow 
velocity is probably best done from a combined analysis of our interferometry measurements 
and single particle spectra. 

At central rapidity the K~/K + ratio falls as the baryon density increases because of 
the increasing scarcity of the anti-up quarks required to form K~ . The increasing K + /tt + 
ratio as well as the constancy of ^ ' K+ p .p implies that the fraction of strange quarks in the 
hadronic system is increasing. 
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TABLES 



TABLE I. Fraction of events satisfying the interaction trigger, target thickness and number of 
events at each setting for each particle type. The target thickness, A, is quoted in nuclear collision 
lengths for the given system. 





Event 


A.ngle 


\ 


4GeV/c 


8GeV/c 


4GeV/c 






Fraction 


mrad 


% 


K+ 


K~ 


K+ 


K~ 


7T+ 


7T+ | 


pBe 


84 ± 1.5% 


44 


3.3 


22723 


47498 


3104 


11594 


15583 


69088 






131 


3.3 


5801 


7948 


2673 


2671 


38240 


15677 


pS 


90 ± 2 % 


44 


3.3 


40389 


39161 


1657 


11263 


15594 


67043 






131 


3.3 


22949 


4661 














pPb 


97 ± 3 % 


44 


4.7 


9591 


45303 


1304 


12268 


13147 


22686 






131 


9.9 


8898 


1304 


4674 


3497 


25553 


5195 


SS 


8.7± .5% 


44 


6.6 


9475 


11875 


6539 


10236 


20201 


21155 






131 


6.6 


9907 


9996 


2780 





25898 


16777 


SPb 


10.7L .6% 


44 


5.9 


8893 


8695 





14034 


3668 


4613 






131 


5.9 


18653 


20028 


3645 





22667 


15949 



TABLE II. Systematic errors on the inverse slopes and dN/dy. 







Kaons 






Pions 




System 


Error on 


Inverse Slope 


Error 


Error on 


Inverse Slope 


Error 




y=2.3-2.9 


y=2.4-3.5 


on dN/dy 


y=2.4-3.0 


y=3.1-4.0 


on dN/dy 


pBe 


10 MeV/c 


10 MeV/c 


9% 


10 MeV/c 


10 MeV/c 


9% 


pS 


10 MeV/c 


10 MeV/c 


9% 


10 MeV/c 


10 MeV/c 


9% 


pPb 


10 MeV/c 


10 MeV/c 


10% 


10 MeV/c 


10 MeV/c 


10% 


SS 


20 MeV/c 


10 MeV/c 


9% 


10 MeV/c 


10 MeV/c 


8% 


SPb 


20 MeV/c 


10 MeV/c 


14% 


10 MeV/c 


10 MeV/c 


9% 
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TABLE III. Inverse slopes (in MeV/c) for K + , K and ir + with statistical errors. 





y=2.3-2.9 


y=2.4-3.5 


y=3.1-4.0 


y=2.4-3.0 




k T < 0.33 GeV/c 


k T < 0.84 GeV/c 


k T < 0.64 GeV/c 


k T = 0.3 - 1.2 GeV/c 




K+ 


K~ 


K+ 


K~ 


7T+ 


vr+ 


pBe 


138 ±7 


145 ±8 


154 ±8 


153 ±4 


148 ±3 


169 ±3 


pS 


152 ±4 


154 ± 10 


163 ± 14 


160 ± 12 


139 ±3 




pPb 


151 ±6 


147 ± 18 


172 ±9 


152 ±5 


145 ±3 


167 ±2 


SS 


181 ± 7 


177 ±6 


159 ±5 


163 ±5 


154 ±5 


179 ±4 


SPb 


197 ±13 


213 ±10 


181 ±8 


175 ±4 


165 ±9 


192 ±5 



TABLE IV. dN/dy for K + , K and ir + with statistical and systematic errors. 





y=2.3-2.9 


y=2.4-3.5 


y=2. 4-4.0 




K+ 


K~ 


K+ 


K~ 




pBe 


.18 ± .01 ± .02 


.14 ±.01 ± .01 


.15 ± .02 ± .02 


.13 ± .01 ± .01 


2.1 ±0.2 ±0.2 


pS 


.24 ± .01 ± .02 


.21 ± .01 ± .02 


.24 ± .03 ± .02 


.18 ±.02 ±.02 


2.4 ±0.2 ±0.2 


pPb 


.29 ± .01 ± .03 


.22 ± .02 ± .02 


.27 ± .02 ± .03 


.21 ± .02 ± .02 


3.2 ±0.2 ±0.3 


SS 


3.76 ± .19 ± .35 


2.4±.1±.2 


4.9 ± .4± .5 


3.4 ±.3 ±.3 


29.7 ± 1.8 ±2.2 


SPb 


9.0 ±.5 ±1.2 


5.4 ±.2 ±.7 


13.1 ± 1.6 ± 1.8 


6.7 ± .4± .9 


62 ± 7 ± 5 
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Cherenkov Beam Counter 
Beam Veto Counter 



FIGURES 

TO Counter 

Silicon Multiplicity Counter 




Beam direction 



MWPC3 




MWPC4 



UCAL 



D: Dipole Magnet 



C: Threshold Gas Cherenkov Counter 
H: Scintillator Hodoscope 
MWPC: Multi-Wire-Propotional Chamber 



Q: Superconducting Quadrupole Magnet UCAL: Uranium-Cu-Scintillator Calorimeter 

FIG. 1. The NA44 experimental set-up. 
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FIG. 2. The acceptance of charged pions and kaons in y and px- The 
decreases from 2tt at pt = 0.0 GeV/c to 0.1 at pt = 0.8 GeV/c. 
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FIG. 3. Mass-squared distribution from Hodoscope 3 with Cerenkov 1 in coincidence for 
(top) and veto mode for kaons (bottom). 
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FIG. 8. The inverse slopes of pions kaons and protons vs system, the errors are statistical. 
Correcting for weak decays would increase the proton inverse slopes by 2 ± 1% for pBe and 10 ± 2% 
for SPb, see [11]. 
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FIG. 10. The K /K + ratio vs system for data and RQMD. Systematic and statistical errors 
have been added in quadrature. 
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FIG. 12. K /K + versus p/p for various systems at y = 2.3 — 2.9. 
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FIG. 13. The K + jp and K jp ratios versus system. 
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